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■^i- Abstract. We establish an effective Markov theory for the rotational Brownian motion of 

y—i hot nanobeads and nanorods. Compact analytical expressions for the effective temperature 

and friction are derived from the fluctuating hydrodynamic equations of motion. They are 
ff~ verified by comparison with recent measurements and with GPU powered parallel molecular 

dynamics simulations over a wide temperature range. This provides unique insights into the 
physics of hot Brownian motion and an excellent starting point for further experimental tests 
C$ and applications involving laser-heated nanobeads, nanorods and Janus particles. 
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1. Introduction 

The popular Markovian theory of Brownian motion, as developed by Einstein, Langevin, 
and Smoluchowski a century ago, has been the starting point and inspiration for innumerable 
applications [Hl|2||. However, the usual convenient formulation in terms of the centre-of-mass 
coordinates of particles only pertains to the special case of an isolated spherical particle. In 
the general case of interacting and/or anisotropic particles, both translational and rotational 
degrees of freedom couple, calling for a more elaborate mathematical description. This is 
most obvious for rod-shaped particles that have different mobilities for the movement parallel 
and perpendicular to their long axis [0, but in fact also holds for interacting spherical particles 
PI . Due to the associated technical complications, the present theoretical understanding is 
still relatively incomplete [5 1, in particular with regard to micro- swimmers and other active or 
self-propelled colloidal particles (6l Ul ID, f° r which the proper hydrodynamic description is 
even more subtle than for passive particles in external fields D9l[T0l. The directed motion for 
such self-propelled particles from sperms |[TT| to Janus particles running on chemical fuel |fT2j 
is usually limited by (equilibrium or nonequilibrium) rotational Brownian motion. Besides, 
rotational Brownian motion is undoubtedly of interest for its own sake. It is accessible 
to spectroscopy |fl3j and has been the basis for the development of new microrheological 
techniques IPT41 and nanoscopic heat engines [ 15 1. 

In this paper, we are concerned with a specific type of rotational Brownian motion that 
occurs whenever the colloidal particles have an elevated temperature with respect to their 
solvent. In this case, we speak of rotational hot Brownian motion, in analogy to the better 
understood translational case ifTBll . Both intended ifTTl [T8l [T9l l20l |2TT| and unintentional 
Il22l |23l realizations of (rotational) hot Brownian motion are nowadays widespread in 
biophysical and nanotechnological applications, which often employ nanoparticles exposed to 
laser light as tracers, anchors, and localized heat sources. Deliberate heating of nanoparticles 
is, for instance, common in photothermal therapy fl24l 1251 . but it also helps to enhance 
the optical contrast for detection [26] or in photothermal correlation spectroscopy [27J. 
Laser-heating is also a convenient way of supplying the energy for the self-thermophoretic 
propulsion of anisotropic particles [fT7|. A quantitative theory for the optical scattering from 
dissolved hot nanoparticles has only recently become available 11281 , paving the way for a 
broad range of future applications. 

In the following, we show that the heating affects the rotational and translational 
degrees of freedom differently, which is due to the intrinsic nonequilibrium nature 
of the phenomenon. To this end, we perform analytical calculations based on 
nonequilibrium fluctuating hydrodynamics, which we compare to large-scale molecular 
dynamics (MD) simulations of an atomistic model of a nanoparticle dissolved in a 
Lennard- Jones fluid. Although, in practice, even for an isolated single colloid rotational 
and translational diffusion always occur simultaneously, we may focus on the one or 
the other separately, in theory. The requirement for this considerable simplification 
is that the coupling only results in a superposition of the respective displacements in 
space and orientation. To be more specific, we assume that the conditions governing 
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the rotational Brownian motion, i. e., the spatially heterogeneous solvent viscosity and 
temperature around the particle, do not depend on the translational Brownian motion. 
This assumption relies on the common Brown- 
ian scale separation. Typical nanoparticle diffu- 
sivities are on the order of 10~ u . . . 10~ 10 m 2 s _1 
while heat and vorticity diffuse at 10~ 7 m s -1 . 
Thus, the hot Brownian motion of a single spher- 
ical nanobead constitutes a stationary nonequilib- 
rium process with fixed radial temperature and 
viscosity fields, T(r) and Tf(r), respectively, if 
the instantaneous particle position is taken as the 
origin of the coordinate system. On this basis, 
we construct and validate a Markov model for the 
rotational Langevin dynamics of a hot Brownian 
particle with effective temperature and friction 
parameters, T^ BM and £hbm- While the success 
of this strategy has already been demonstrated 
for the translational motion lfT6l l29l l30ll . recent 
experiments using heated nanorods ||23l and hot 
Janus particles ifTTl underscore the need for a sep- 
arate quantitative analysis of the rotational dynamics. Below, we derive T^ iBM and Chbm f° r 
the rotational Brownian motion of a hot particle and demonstrate that they differ from their 
analogues for translational motion. The mathematical structure of the rotational dynamics is 
simpler and allows for analytical solutions where one has to resort to numerical methods in 
the translational case. 

The effective Langevin equation for the rotational dynamics of the nanoparticle 
orientation n reads PT1. 




Figure 1. Artist's conception of a hot Brownian 
particle, illustrating some notation. 



CHBMn = ^ xn - (!) 

Under the presupposed nonequilibrium steady-state conditions outlined above, £hbm ls me 
effective rotational friction coefficient. In the isothermal limit, £hbm _ > 87TT]7? 3 for a sphere 
of radius R in a solvent of viscosity r\ . The stochastic torque £> is assumed to be a Gaussian 
random variable characterized by the moments 

(5(0) =0, <6(0&(O> = 2fe7S HM 8 BM fy5(f-0 . (2) 

The form of the noise strength amounts to the assumption that a generalized Einstein relation 

^bThBM = ^HBM bHBM (3) 

links the effective friction and temperature, Chbm an( ^ ^HBM' t0 tne effective rotational 
diffusivity £>hbm- That the Brownian dynamics of a single heated particle is indeed 
constrained by this quasi-equilibrium relation is corroborated by our MD simulations, 
presented below. 



Rotational hot Brownian motion 4 

2. Theory 

The classical rotational Stokes problem is to find the friction coefficient of a steadily rotating 
particle in a viscous fluid of homogeneous solvent viscosity. For heated particles, the 
assumption of a constant viscosity has to be relaxed. The general case of an arbitrarily shaped 
hot particle, which induces an asymmetric temperature profile in the solvent, gives rise to 
formidable technical complications. With the aim of deriving analytical results, we restrict 
our discussion to spherical beads (cf. figure [T]). Then the temperature can be idealized as a 
radial field T(r) that entails a radially varying viscosity r\ (r) via some constitutive law of the 
solvent, which we assume to be given. Some extensions to spheroids and slender cylinders 
can be discussed along the same lines if the radius r is given a slightly different interpretation, 
as outlined below 11321 . We further take the solvent to be incompressible, which is a good 
approximation for most common solvents, such as water, and eases the calculation. In the 
low Reynolds number limit, applicable to micro- and nanoparticles in solution, the solvent 
velocity u(r) follows from 

vp = 2v-T]r = 2rvT]+2T]V-r, v- u = o, (4) 

with the strain rate tensor V = (Vu + Vu r )/2 and the pressure p. 

For the purely rotational fluid motion around a steadily rotating sphere, analytical 
solutions can be found, as follows. The high symmetry of the velocity field 

u(r, 8,<j>) =w^e^ =w(r)sin0e0 (5) 

around the sphere entails a highly degenerate strain rate tensor. The tensor element Tq^ 
contains u§ only in the form d$ (us/ sin 6) = dgw(r) = 0, leaving us with only a single relevant 
entry T^, r with 

2T§ r = d r U(j, — Uty/r = (w — w/r)sin0 . (6) 

Since the viscosity is assumed to vary only radially, its gradient is V?7(r) = rj'e r , and the 
equation of motion, equation (|4]), reduces to 

Vp = [2t]TV + tj V 2 K - i\u$r~ 2 sin" 2 0] e (7) 

Due to the cylindrical symmetry, the pressure gradient must not contain an azimuthal 
component, though, i. e., Vp = 0. This leads to the ordinary differential equation for w(r), 

w" + (2/r + ri'/ri)(w'-w/r)=0. (8) 

For the known case of a constant viscosity 7] (r) = 7]q one easily verifies by insertion the 
solution w = c\r-\-C2/r 2 with constants c\ and c%. For an unbounded fluid c\ = so that the 
flow field takes the familiar form u = Ox r(R/r) 3 for a sphere of radius R rotating at constant 
angular velocity Q. A constant viscosity can be interpreted as a degenerate case (for n = 2) of 
the power law viscosity field 

r,(r) = T 1 o(r/R) n - 2 . (9) 
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For the latter T\' jt\ = (n — 2)/r, hence equation ^ reduces to w" + nw'/r — nw/r 2 = 0, which 
has the general solution w(r) = c\r + C2r~ n . Thus the flow field in an unbounded fluid with a 
power-law viscosity field with arbitrary n > 0, reads 

u = Qxr(i?/r)" +1 . (10) 

For more complex viscosity profiles 7](r), the task of solving equation (J8J) can be reduced to 
an integration using the Wronskian. Knowing the particular solution w(r) = r, the general 
solution to equation ([8]) is found to be 

W (r)= Cl r + c 2 rJ t ^i-,dx. (11) 



Moreover, rewriting equation d8j) in terms of y = w' — w/r yields 

Y + ([lnr 1 ]' + 3/r)y=0, (12) 

which is (up to a constant factor) solved by 

y(r)oc(/?/ r ) 3 /T?(r). (13) 

With the solution of the Stokes problem at hand, we can explicitly calculate the effective 
friction coefficient Chbm f° r a rotating sphere. The torque exerted by a rotating sphere on the 
surrounding fluid is obtained by integrating IrfT^r-r xe^ = rjys'mOr x e* over the surface 
of the sphere, where r x e^ = — R&q = Rs'm6e z — Rcos6e p . As the (^-integral over the e p - 
component vanishes identically, this yields 

fit Qjt 

In d6 R 3 sin 3 6 tj (R)\y(R)\ = —R 3 tj (R)\y(R)\ , (14) 

Jo 3 

and division by \Cl\ yields the wanted effective friction coefficient. The boundary conditions 

u(r)|| r |=R and limi r i^ 00 u(r), and the viscosity profile r\{r) fix the constants c\ and c 2 in 



equation (11). For an unbounded medium, w(r — > °°) = implies c\ = 0, and w(r = R) = CLR 



fixes C2- Collecting results, we find 

O t-ca 1 

«" =rJ R W)S ir ' <15) 

which indeed attains the isothermal value l/(87rr/o^ 3 ) for r\ (r) = t\q. 

Following the derivations in Ref. [30J, the second important parameter of a hot Brownian 
particle, its effective Brownian temperature ?hbm> is given by 

J y r(r)0(r)d 3 r 
r ™ M= J^(r)d 3 r • (16) 

(Depending on the type of motion of the particle, the result gives the translational/rotational 
effective temperature T^ M .) It can explicitly be determined by integration once the 
dissipation function = rj V : l~/2 is known. For the rotational flow field in equation ([5]), 

<j) = ri(w' -w/r) 2 sin 2 = Tjy 2 sin 2 6, (17) 

from which we get 

HBM " £ri-i(r)r-4dr ' (18) 
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Together with equation ( 15 ), this completes our formal derivation of the effective friction and 
temperature parameters characterizing the rotational hot Brownian motion of a sphere. They 
can explicitly be evaluated and used in the effective Langevin equation ([T]) & ([2]), provided 
that the temperature dependence 7] (T) of the solvent viscosity and the temperature profile 
T(r) around the nanoparticle are known. For the practically important special case that the 
solvent is water, which is well characterized by 

t\{T) = r] 00 e A/{T - TYF) and T(r) = To + ATR/r, (19) 

the convenient approximation (accurate to within 2% for AT < To) 

Thbm « To + 3 -AT (20) 

is obtained by neglecting the temperature dependence of the viscosity, i. e., by setting 



r\ (r) = constant, in equation ( 18 ) 



For the example of a Lennard- Jones fluid with an immersed hot nanobead, we have 



evaluated equation ( 18 ) numerically, based on the temperature and viscosity profiles obtained 
from our simulations. Care has been taken to include the density variations and finite size 
effects correctly, as detailed in Ref . [J30!] . In figure 15} we compare the predicted T^g M to the 
effective temperature deduced from the directly measured friction Chbm and diffusively D I HBM 
via the generalized Einstein relation Q. (A brief description of the MD simulations can be 
found further below.) 

Beyond the above limiting results for spherical particles, analytical estimates for the 
rotational hot Brownian motion of anisotropic particles can readily be obtained in the slender 
rod limit. To find the friction coefficient per unit length for an infinitely long hot cylinder, 
equation Q has to be extended by adding Oseen's term pU- Vu to the force on the right hand 
side [33] . The effective friction coefficient for stationary translation along the main axis can 
then be calculated as [34J 

£1bm = 2k| / - »,L dr| , (21) 



e 

r\{r)r 



where the radius r* now denotes the distance to the symmetry axis of the cylinder. Exploiting 
the formal analogy between the Reynolds number pUR/rj (R) and the aspect ratio L/(2R) in 
cutting off the hydrodynamic divergences under isothermal conditions, we arrive at a plausible 
estimate for the friction coefficient of a slender cylinder of finite length L and radius R, namely, 

4 BM = 27tr,(R){\n[r,(R)L/(2r, R)] -y)" 1 . (22) 

For transverse and rotational motion, expressions corresponding to equation pT| ) are more 
complicated to calculate, due to the lower symmetry of the temperature and fluid velocity 
fields. However, simple estimates for the effective transverse and rotational friction 
coefficients are readily obtained if one neglects a potential dependence of the ratios of the 
various friction coefficients on the heating, namely, 

Chbm ~ 2 Chbm and "=hbm ~ ( L / 12 )?hbm ■ ( 23 ) 
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Figure 2. Left: Effective temperature of rotational hot Brownian motion. The simulation 
results for T^ BM (•) were deduced from the numerically measured Chbm anc ' ^hbm usm g me 
generalized Einstein relation Oil. An alternative estimation of 7j"i BM (□) was obtained from 
the Boltzmann distribution of the inclination angle 9 in a harmonic angular confinement (right 
panel). The theoretical prediction (solid line) was evaluated within the idealized theory for an 



incompressible fluid via equation ( 18 1 using the radial viscosity profile 77 (r) determined in the 
MD simulation. For comparison, the effective temperature 7^ BM for the translational degrees 
of freedom (dot-dashed line) [30 1 and the solvent temperature at the particle surface (dotted 
line) are shown. Right: The measured distribution of the inclination angle 6 in a harmonic 
angular confinement potential for nanoparticle temperatures T p = 0.8£/&b (•), 1.25 £/£b (■) 
and the corresponding distribution p(9) ~ <?~P y ( e ) with /3 _1 = ^b^hem depicted by the solid 
lines. 



The result for Chbm' normalized to its isothermal limit, is depicted in the left panel of 
figure [3} for various aspect ratios. For aspect ratios that are large enough to admit the 



slender-rod approximation, the temperature dependence of the expression (22), normalized 
to its isothermal limit, is close to the normalized rotational friction coefficient of a sphere. 

A practical estimate for T^ BM of a hot rotating rod is obtained to first order in AT 
by ignoring the temperature dependence of the viscosity and treating the rod as a prolate 
spheroid. In this case, the flow field u(r) is known analytically [35] . It is constructed by a 
line distribution of some moments of the fundamental solutions to the Stokes equations under 
different singular forcings, so-called stresslets, rotlets, and potential quadrupoles. From u(r), 
a straightforward calculation yields the "radial" dissipation function (t) from which T^ BM is 
calculated for a given temperature field T(r) around the spheroid. In spheroidal coordinates, 
(x, C,<P), where x = cxC,, y = ca/(t 2 - 1)(1 - £ 2 )cos<p, z = ca/(t 2 - 1)(1 -£ 2 )sin<p, the 
temperature field 

T ( t) = Tq + ATarccot r/arccot t (24) 

obtained from Fourier's law for constant heat conductivity only depends on the "radius" x and 
the eccentricity Xq 2 = 1 — (2R/L) 2 of the particle. Using it together with the isothermal 0(t) 



in equation ( 16), we obtain the first order term of the series expansion of T^ BM for slender 
particles in powers of AT. Its coefficient (T^ BM — Tq) j AT is plotted in figure [3} Deviations 
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from the value 3/4 for a sphere, equation ( 20 ), only become apparent for extreme aspect ratios 
L/(2R) > 10 2 , corresponding to eccentricities Tq 2 > 0.9999. 

Note that the high power of the radial distance r in the denominator under the integral 



in eqs. (15) & (18) suggests that the effective rotational temperature and friction should 
be higher and lower than their translational counterparts, respectively. This conclusion is 
supported by recently published experimental data for gold nanorods, where the effective 
temperature for their rotational Brownian motion in an optical trap is "found to be close to the 
particle's temperature" [|23l . The plots of our predictions in figure [3] suggest that our results 
for a spherical particle should still provide a reasonably good approximation for these rods, 
which have an aspect ratio of about two. To make closer contact with the experiments, we 



have used equation ( 15 ) and the differential shell method [12911 to calculate the rotational and 
translational friction coefficient for the viscosity field of equation ( 19). The inset of figure [3] 
depicts the ratio of the estimated rotational and translational relaxation times T rt <* Q BM /AT 
for the experimental particle in the optical trap, assuming a linear dependence between the 
trap stiffness and the heating AT. The prediction is seen to be in good qualitative agreement 
with Fig. 3 of Ref. 11251. 

One may ask, whether any more general statements can be made concerning the relative 
magnitude of the various effective temperatures of hot Brownian motion introduced so far, 
independent of the shape of the particle, and based on some generic material properties (e.g. 
that the solvent viscosity decreases upon heating). In particular, one expects that, under 
otherwise identical conditions, the effective temperature T^ BM will usually be higher for 
rotational than for translational motion, because the solvent velocity field around the particle 
is more localized near the hot particle for rotation than for translation. This suggests that the 
temperature ordering 

To < Trbm < ^hbm <T S <T V (25) 

might be quite generic (T p and T s are the particle temperature and the solvent temperature at 
the particle surface, respectively, Tq is the ambient temperature). The claim is corroborated 
by the results for a temperature-independent viscosity and for a viscosity step (see appendix), 
and also by experimental observations lTT6ll23l . 

3. MD simulations 

Numerical simulations allow for a more accurate check of some of our theoretical prediction 
than the quoted experiments, since we can better control the "experimental" conditions. 
Our simulations of rotational hot Brownian motion are based on the same setup as in the 
translational case 11301 . In brief, the system is modelled as a Lennard- Jones fluid with a 
radial pair potential V(r) = 4e[(a/r) 12 — (o/r) 6 ]. While 107233 particles are treated as 
solvent, a spherical cluster of 767 particles, additionally interconnected by a FENE potential 
V(r) = -0.5fd?olog[l - ( r / R o) 2 } w i tn K = 30e/a 2 , Rq = \.5a, forms the nanoparticle. Each 
simulation run consists of an initial isothermal equilibration phase at T = 0.75 e/k^, followed 
by a heating phase, where the nanoparticle is kept at a constant elevated temperature T p , and 
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Figure 3. Left: Effective friction coefficient £hbm of a hot nanoparticle in water (normalized 
to the isothermal limiting value) as a function of the temperature increment Ar above the 
ambient temperature To = 298 K. The analytical predictions for a spherical bead of radius R 



from equation ( 15 1 and Ref. |29| are represented by solid and dot-dashed lines for rotational 



and translational motion, respectively. Note that they differ by a heating-dependent kinematic 
factor, an effect that we neglect in the estimate equation ( |23] >, for slender rods. Dashed lines 
represent the effective longitudinal friction coefficient Chbm °f a ^ ot s l en der rod according to 
equation (22 1, for the aspect ratios L/(2R) = 10 1 , 10 2 , 10 3 , 10 4 . The curves closely follow the 
prediction for the rotating sphere, except for small L/(2R) < 10 2 near the isotropic limit, which 
is not correctly recovered by the slender-rod approximation. Inset: The inverse translational 
versus the inverse rotational relaxation time of a spherical particle in an optical trap of 
varying strength; to be compared with Fig. 3 d) of Ref. IT231 . The dotted line indicates the 
naive estimate obtained by identifying the effective viscosities for translational and rotational 
motion. Right: Variation of the first order coefficient of T^ BM (AT) with the eccentricity 
Tg of the spheroidal particle. No numerically stable integration was attained that covers the 
whole range 0. . . 1. Therefore, a 7th order series expansion of the dissipation function 0(t) 



-2 



in T was employed to generate the correct asymptotic behaviour as zZ — > 0, while the full 
expression for 0(t) was evaluated numerically for the other branch. As a guide to the eye, a 
matching curve has been added by hand. 



a boundary layer at the periphery of the simulation volume is kept at T = 0.75 £/k&. All 
measurements were performed once a stationary state had been attained in the heating phase. 
The results were averaged over 30 independent trajectories of 2 x 10 5 steps (corresponding to 
a physical duration of Ins). For further details concerning the simulation method, we refer 
the reader to Ref. W& 

The rotational diffusion coefficient D r WM was determined from the decay of the 
autocorrelation function of an orientation vector assigned to the freely diffusing hot 
nanoparticle figure [T| The rotational friction coefficient £hbm was determined by applying a 
constant external torque to the nanoparticle and measuring its angular velocity. To gain more 
insight into the microscopic details, we also recorded the radially varying angular velocity 
in concentric shells of thickness a/5 around the nanoparticle, in this case. Figure [4] shows 
the excellent agreement of the measured angular velocity field w(r)/r with the theoretical 
prediction in equation (111- except in the close vicinity of the nanoparticle (inset of figure [4]) 
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Figure 4. Radial variation of the measured angular velocity w(r)/r normalized by the angular 
velocity of the hot nanoparticle for T ? = 0.9e/k^ (•) and 1 .25e/A:b (')• The solid lines are the 
corresponding plots of equation ( |TT[ > evaluated using the measured viscosity and temperature 
profiles. Inset: The variation of the angular velocity near the nanoparticle surface. The vertical 
solid lines indicate the positions of the hydrodynamic boundary condition for the predicted 



hydrodynamic flow fields (Hi fitted to the simulation data. The shaded region (its extension 
marked by an arrow) indicates the equivalent sphere radius for the nanoparticle (dot-dashed 
line) plus a/2. The comparison reveals a weak apparent slip at high nanoparticle temperatures 
due to the radial solvent density variation induced by the heating. 



where solvent density variations are prominent, which have been neglected in the analytical 
calculation. For a comparison with the predictions, we therefore treated the radius R of the 
nanoparticle, which affects w(r) via the boundary condition, as a free parameter. When the 
heating of the nanoparticle is small, R agrees with the value Ru + a/2, with the radius 7?h of 
an equivalent solid sphere given by R^ = gjL £f ; - =1 (r,- — r ; ) 2 [361 137]|. For strong heating, 
R decreases slightly, due to the mentioned solvent density variation, which results in a weak 
surface slip of the hydrodynamic flow field (inset of figure [4]). 

Finally, we pursued three routes for determining the hot nanoparticle's effective 
temperature: (i) by measuring the rotational diffusion and friction coefficients, D r nBM and 
ChbM' an d making use of the generalized Einstein relation in equation ^ (• in figure |2l); 
(ii) by fitting the recorded angular distribution p(6) in a harmonic angular confinement 
potential V{9) = f (0 - n/2) 2 with a Boltzmann distribution p(6) °c exp(-V(0)/fc B 7^ BM ) 
(□ in figure [2]); and (Hi) by measuring the viscosity and temperature profiles, r](r) and T(r) 
and calculating 7^ BM via equation (18) (solid line in figure |2l). All three methods agree 
within the error bars, thereby confirming our theoretical predictions for the effective friction 
£hbM' me effective temperature Tj^Bivi' an( ^ tne generalized Einstein relation of rotational hot 
Brownian motion, equation ([3]). 



Rotational hot Brownian motion 1 1 

4. Conclusions 

We have analyzed the rotational Brownian motion of a heated (nano-)particle in an otherwise 
unheated solvent. Particles of spherical, spheroidal and slender cylindrical shape have been 
considered. While we have obtained exact results for the effective rotational friction and 
temperature in the case of spherical particles, we contended ourselves with approximate 
estimates for prolate shapes. Disregarding a potential temperature dependence of the ratio 



between the transverse and the longitudinal friction for a rod in equation (23) is justified for 
large particle aspect ratios, since the near-field solutions for the solvent flow are dominated by 
terms with the same radial dependence. Also, while Oseen's treatment of the infinite cylinder 



in equation (21 ) adds a nonlinearity in U to the Stokes problem, calling for a nonlinear (and 
non-Markovian) theory, the linear analysis of Ref. [30] remains valid for finite particles. It 
leaves only the (possibly difficult) technical problem of evaluating the flow velocity field from 
which the effective parameters are calculated. 

We validated our theoretical predictions against large-scale molecular dynamics 
simulations of a Lennard- Jones system with an immersed hot colloid, and by comparison 
with recent measurements of diffusing hot nanorods [|23l . Thereby, we could quantitatively 
establish an effective equilibrium description for the nonequilibrium rotational motion of 
a heated Brownian particle. While the general procedures and results resemble those for 
translational motion 11291 . their analytical tractability and their detailed structure was found 
to be substantially different. From a comparison of the rotational and translational case, an 
ordering of the effective temperatures governing hot Brownian motion emerged. Our findings 
provide the basis for a quantitative investigation of the Brownian dynamics of particles of 
anisotropic shapes, interactions between hot particles [38|, and self-thermophoretic Janus 
particles that move on persistent Brownian paths determined by a mutual competition of 
rotational and translational hot Brownian motion [17J. 
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Appendix: Stokes' problem for a viscosity step. 

From equations (49) and (53) of Ref. [|29l . the viscous dissipation function 0(r) = 
/ 0(r) dOdcp for a translating sphere, driven by the constant force F = Fe z , follows as 

0(r) = ^ [l5al-3Fa 3 r 2 /(2nri) 

+ r 4 (F 2 /(16n 2 T] 2 )+F/(nri)a2r 3 + 10a 2 2 r 6 ^ , (A.l) 

which depends on r\(r) in an intricate way, namely via equations (44) of Ref. [29] . While 
this allows for an explicit computation of 7h BM once r\(r) is specified, general statements 
about the relative magnitudes of 7^ BM and T^ BM seem thus hard to attain, even for the 
highly symmetric case of a spherical bead. We therefore contend ourselves with outlining the 
principles how such general considerations might proceed, namely by considering a spherical 
bead and the class of step profiles, 

t < PR : Tj(r) = Tjo/K , r > PR : Tj(r) = tj , (A.2) 

for which both the translational and rotational effective temperatures can be calculated 
analytically. Note that the common case of a fluid viscosity that decreases upon heating 
corresponds to K> 1, while K < 1 holds for dilute gases. 
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The effective rotational temperature 7^ BM follows from equation (18) and takes the 
simple form 

3(/3 4 -1)k + 3 



[ HBM 



To 



ATR 4(/3 4 k-/3k + j8)' 

where Tq is the ambient temperature. 

For the translational case, the general ansatz for the velocity and pressure fields, u 
u r e r + ueeo and p, reads [39] : 



(A.3) 



a\ 2 ay. 

ao~\ Ya 2 r + — Icost* 



«3 



Ur{r,6) 

u e (r,0)- ru , ^ , wz . ^ 

p(r, 6) =p +(-^ + 10a 2 r)ricos6 






sin0 



(A.4) 
(A.5) 
(A.6) 



where u r and we denote the radial and polar velocity components, respectively. While the 
coefficients are l l39l 

(A.7) 



ao = uq, a\ = —3uqR/2, a 2 =0, a 3 = uqR /2 



in an infinite homogeneous system, where uq denotes the particle velocity relative to the 
resting fluid at infinity, the coefficients at each take two r-dependent values for the step profile 
of the viscosity. These coefficients are found by matching the boundary conditions for the 
velocity and the stress at r = R, r = J5R, and r — > oo. More precisely, the following linear 
system of equations needs to be solved: 



with 
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Ma 
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-2R 3 

R 3 p 3 
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(A.8) 
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(A.9) 



a 2 (r> pR)=0. 



(A.10) 
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We find as solution, 

/ /3k[/3 5 (4k + 6) + 5/3 2 (k-1)-9k + 9] \ 

-3R[5k[j5 5 (2k + 3)-2k + 2] 

-3pKR- 2 (p 2 -l){K-l) 

;? 3 /3 3 k[/3 3 (2k + 3)-2k + 2] 

-3Rp[p 5 (2K + 3)-2K + 2] 
\ 7? 3 /3 3 [-3/3 5 (k-1)+5/3 3 k-2k + 2] / 

where <€ = 4 + 6/3 5 + (6/3 6 + 3/3 5 - 10/3 3 + 9/3 - 8) K+ (j8 - l) 4 [jS (4/3 +7) +4] k 2 . 



[7 



(A.ll) 



We use this solution together with equation (A.l ) to compute 



J HBM 



rj5R 



fg ^rdr+[g R ^rdr 



ATR 



rpR 



/jf (j)r 2 dr + Jp R ^r 2 dr 



(A.12) 



Dividing the expression in equation ( A.3[ ) by this result and evaluating the integrals, we obtain 

T ° '""• - "-- "•'" (A.13) 



J HBM 
Tt 



_ _ PQ+PlK + p 2 K 2 + p 3 K 3 +p 4 K 4 
Tq~ qo + qiK + q 2 K 2 + q 3 K 3 + q 4 K 4 



where (/3 = 1 + e) 

Po-qo 

pi-qi 



s(3(e + \) 5 +2^ 

e(l62e 13 + 2178e 12 + 13482c 11 +51030e 10 + 131931c 9 
+ 246390c 8 + 343605c 7 + 366360c 6 + 304685c 5 



+ 200610c 4 + 104325c 3 + 41350c 2 + 1 1400c + 1800) 
Pi- qi = 3e 2 (24e 13 +423e 12 + 3252c 11 + 14743c 10 + 44630e 9 

+ 96664c 8 + 156760c 7 + 196480c 6 + 193820c 5 

+ 150415c 4 + 89750c 3 + 39325c 2 + 1 1550c + 1800) 
p 3 - q3 = e 3 (96e 12 + 1251c 11 + 7704c 10 + 30021c 9 + 83280c 8 

+ 174165c 7 + 281090c 6 + 349995c 5 + 330540c 4 

+ 229950c 3 + 1 1 2500c 2 + 35 1 00c + 5400) 
p A - q A = 2e 8 (l6e 7 + 159e 6 + 681e 5 + 1634c 4 + 2400c 3 + 2220c 2 

+ 1260c + 360) . 

Obviously, since /3 > 1 by definition, and thus e > 0, each coefficient p\ in the numerator 
of equation ( |A.13[ ) is strictly larger than the corresponding coefficient in the denominator, g,-. 
Hence, for a viscosity profile that increases from r\ (R) to T]o in a single step, we have shown 
that 7gu M > Z£jj M holds. 
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